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^ , Abstract. 

The discrete equations of motion for the quantum mappings of KdV type are given in terms of 
00 ■ the Sklyanin variables (which are also known as quantum separated variables). Both temporal 

(discrete-time) evolutions and spatial (along the lattice at a constant time-level) evolutions are 
considered. In the classical limit, the temporal equations reduce to the (classical) discrete 
^ _. , Dubrovin equations as given in a previous publication (Nijhoff F W 2000 Chaos, Solitons and 

• ' Fractals 11 19-28). The reconstruction of the original dynamical variables in terms of the Sklyanin 

" ■ variables is also achieved. 
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'nI" ' 1. Introduction 

O 

^ . 

O ■ The quantization of discrete-time systems is an outstanding problem within the wider context of 

Q ! the quantum theory of dynamical systems. From the perspective of quantum chaos theory 

(quantum chaology), quantum mappings have been studied by Berry et al. [H]. Quantum 

mechanical systems evolving in discrete-time have also been considered by Bender et al. |3], [H], 

where the emphasis was on the application of finite-element methods. Both of these approaches 

deal with discrete-time systems that, in the classical limit, generically exhibit chaos. 

Integrable classical mappings, as opposed to mappings that exhibit chaos, have 

.^ \ been systematically constructed and studied in recent years, see, for example, [TT], 

[2n] , |2E] , |S21 , [SB] , ISlj ■ A specific set of examples are the so-called mappings of KdV type [Hj, 

, which were classically constructed from periodic initial value problems on the lattice KdV 

partial difference equation. They are the principal model of interest in this paper. Integrable 
mappings are highly nontrivial; however, there is a great deal of understanding of them, both on 
the level of their solvability as well as on their classification. 

Historically the quantization of integrable models formed many of the paradigms of quantum 
theory. Furthermore, the inherent discreteness of quantum theory points to a specific role that 
discrete, and hence integrable discrete, systems may play in the further development of the 
theory. Quantum integrability on the discrete space-time lattice has been considered in a number 
of papers (mostly from the perspective of i?-matrix theory), for example, [Zj, [Hlj, [20], |2I], |2Z|, 
j35j . Here the discrete aspect is not seen as some kind of approximation but, rather, it is 
postulated as the underlying structure of space-time from the very start. In the present paper. 



Quantum discrete Dubrovin equations 2 

adopting the same point of view, we will be mainly concerned with the approach initiated in [7j, 
[20j . |21j . where a non-ultralocal Yang-Baxter (i?-matrix) structure appropriate for obtaining an 
"integrable quantization" of the mappings of KdV type was given (in the continuous-time setting 
such a non-ultralocal Yang-Baxter structure had been previously used in connection with the 
quantum Toda theory P|). 

Quantizing a discrete-time system is essentially different to the conventional quantization 
procedure of Hamiltonian systems. In the continuous-time setting it is the Hamiltonian 
(or commuting family of Hamiltonians) that stand central to the theory, the spectrum and 
eigenfunctions of which are the main objects to be computed. The equations of motion in terms 
of the canonical operators (in the Heisenberg picture), or the evolution of the states (in the 
Schrodinger picture), only play a subsidiary role. In the discrete-time setting it is no longer the 
Hamiltonian(s) that define the model, but the equations of motion exclusively. This draws us 
away from the conventional schemes of quantization, and leads us to investigate more closely the 
quantum equations of motion. Following the historical imperative once again, integrability offers 
a leading principle to develop our understanding of quantum mechanics, this time in the discrete 
regime. 

By an integrable quantum mapping we mean an automorphism of the quantum algebra 
under consideration, which, furthermore, possesses a "sufficient" set of exact commuting invariant 
functions on the algebra. (We will make this definition more precise later.) This is the quantum 
counterpart of the classical integrability of mappings in the sense of Liouville-Arnol'd-Veselov 
j33j . On the classical level, mappings of this type typically involve rational expressions exhibiting 
singularities that imply that the time-evolution cannot be globally defined. The characterization 
of integrable maps through their singularity structure is a focus point of current investigation, 
see, for example, [I2j, [2H]. In the more restrictive context of the mappings of KdV type, a 
natural resolution is to describe the mapping in terms of variables that live on the Riemann 
surface associated with the underlying spectral curve. The explicit description of the dynamics 
on the Riemann surface was achieved in ^2| (see also [HI, [22]) leading to the so-called discrete 
Dubrovin equations. 

With this insight, one may take the point of view that the proper quantization procedure for 
this discrete-time system is to write the equations of motion in terms of the quantum analogue 
of these variables. These variables are Sklyanin's variables (which are also known as quantum 
separated variables). Following the correspondence with the classical case, the equations of 
motion in terms of the Sklyanin variables are called quantum discrete Dubrovin equations. (It 
should be noted that the Sklyanin variables are taking an increasingly primary role in the field 
of integrable systems. They have played a fundamental part in various recent publications, with 
motivations different to that of this work, such as [l], [2], [3T].) 

The outline of the paper is as follows. The non-ultralocal Yang-Baxter structure of pP is 
recapitulated in section |2l in such a way as to bring the features required for the derivation of the 
quantum discrete Dubrovin equations to the fore. More specific information pertaining to the 
mappings of KdV type is given in sectional In section Hlthe Sklyanin algebra is set up, this paves 
the way for the derivation of the quantum discrete Dubrovin equations in section (Equations 
f)5.14|) and ()5.20|) give a temporal, that is discrete-time, evolution and, hence, are called temporal 
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quantum discrete Dubrovin equations. Equations (J5.29J) and (|5.31|) give a spatial evolution and, 
hence, are called spatial quantum discrete Dubrovin equations.) The reconstruction of the original 
dynamical variables (of the mappings of KdV type) in terms of the Sklyanin variables is also 
addressed in section |S| The well-defined evolution arising from the quantum discrete Dubrovin 
equations and the reconstruction is illustrated in the one and two degrees of freedom situations in 
section ini We remain very formal and algebraic throughout this paper, principally concentrating 
on the derivation of the quantum discrete Dubrovin equations. 

2. Non-ultralocal Yang-Baxter structure 

The non-ultralocal Yang-Baxter structure for the class of discrete-time systems to which the 
mappings of KdV type belong was given in |21j. The convention of that paper (the standard 
convention) will be employed in this section, [Appendix A[ and [Appendix B| This convention 
includes that the subscripts 1,2,... denote factors in a matricial tensor product, and the same 
subscripts distinguish the associated spectral parameters. Care must be taken not to confuse 
these subscripts with those that correspond to the grading of the monodromy matrix or the 
subscripts which identify different dynamical variables. These are, however, all perfectly clear 
within their context. 

The only nontrivial commutation relations between the operators -L„(A) are those on the 
same and nearest-neighbour sites, namely as follows 

Rl2-^n,lLn,2 = -^n,2 -^n,l -^12 ; (2.1a) 

Ln+l,lSi2Ln,2 = Ln,2Ln+l,l , (2-1^) 

Ln,lLm,2 =Lm,2Ln,l \ U - m \> 2 , (2.1c) 

where -L„(A) is the L operator at the nth site and L„j denotes I>„(A) acting nontrivially only on 
the jth factor of the tensor product, 

Lnj := 1 ® 1 O . . . ® L„(Aj) (g) . . . (g) 1. 

jth place 

The operators i?-^ := R-^{\j^ \k) act nontrivially only on the jth and kth factors of the tensor 
product. As was shown in |21j, equations ()2.1aj) to ()2.1c|) lead to 

-"^12 -'n,! -'n,2 — ^ n,2 ^ n,l -"^12 i {2. la) 

^n,l^l2'''n,2 ~ -^ n,2 '-'12 -^ n,l ' [2. 2b) 

where S12 = S21, 



p 



T+(A) := n h^^) T-{X) := n ^.(A) , (2.3) 

j=n+l j=l 
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(P G N denotes the "spatial" periodicity of the model, 1 < n < P — 1) and also lead to 

-^12 ^1 'S'i2 ^2 = T2S^2TlRi2, (2.4) 

where Ti denotes the monodromy matrix, T{X), acting nontrivially only on the first factor of the 
tensor product; the monodromy matrix, 

p 
T(A) := n LnW ■ (2-5) 

n=l 

(The convention for the ordered product in ()2.H|1 and ()2.5j) is such that the matrices L„ are 
ordered from right to left with increasing label n.) In the P = 2 case equation ()2.25|) replaces 
()2.1fe|l . The compatibility relations of equations ()2.1ffl|l to ()2.1c|l lead to the following consistency 
conditions on R"^ and S*^: 

-^12 -^13 -^23 = -^23 -^13 "^12 ) (2.6a) 

-^23 '5^12 '5^13 = '5^13 '5^12 -^23 • (2.6o) 

Equation ()2.6ajl is the quantum Yang-Baxter equation for R^, which is coupled with 5*^ by 
equation ()2.6&|) . It is also assumed that Sj^ and S*^ are invertible. In order to establish that the 
structure given by the above commutation relations allows for suitable commutation relations 
for the monodromy matrix we need to impose in addition to ()2.6flj) and (J2.66J) that 

-^12*^12 = '^12-^12 • (2-7) 

Integrable mappings follow from a discrete-time Zakharov-Shabat system of the form 

Z„(A)M„(A) = M„+i(A)L„(A), (2.8) 

where the tilde, ~, denotes a time update and M„(A) is the discrete-time evolution operator at 
the site n (M„(A) would denote the discrete-time evolution operator at the site n at the next time 
level). The M (or temporal) part of the extended Yang-Baxter structure, as given in [Zj, pOj . 
and |2I] , allows one to derive the invariants of the discrete-time evolution. It also allows one to 
show that the Yang-Baxter relation ()2.4j] is preserved throughout the (discrete) time evolution. 
The only extra relations from the M part of the extended Yang-Baxter structure required for 
these proofs are 

i?+ M„,iM„,2 = M„,2M„,ii?i2 (2.9) 

and 

TiMfi^5+Mi,2 = Mi,2 512 Ti Mf 1^ (2-10) 



The proofs are given in Appendix A 



In the quantum discrete-time setting the commuting family of invariants (i.e., the invariants 
of the discrete-time evolution) are given by expanding 

r(A) = tT{K{X)T{X)) (2.11) 
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in powers of the spectral parameter, A. The proof follows by taking the trace over both spaces 
of the tensor product of P12K1K2 multiplying equation ()2.10j) from the left (where P12 is the 
permutation operator). (The details can be found in [Appendix A[ ) The result is that 

K2 = tT,{P,2''[C^S^2)-']} , (2.12) 

where the left superscript *i denotes the matrix transpose in the first factor of the matricial 
tensor product. Observe that the commutation relation for the L„ operators, equation ()2.1ajl . is 
of the same form as that for the M„ operators, equation ()2.9|) . and it follows immediately from 
equation ()2.2&|) that 

^1 -^1,1 'S'12 -^1,2 = -^1,2 5*12 T^i -^1,1 , (2-13) 

which is of the same form as equation ()2.1()j) . Hence it follows, in an exactly analogous fashion 
to the temporal evolution, that there is a "spatial" evolution which preserves the Yang-Baxter 
relation ()2.4j) and the family of invariants ()2.1H) . This is proven (in both the spatial and temporal 
case) in Appendix A[ The spatial evolution is denoted by the hat,^, hence 

f(A)=Li(A)T(A)Li(A)-\ (2.14) 

For later purposes it is now assumed that R12 is proportional to a rank-one projector for 
a particular relative value of the spectral parameters Ai and A2. This occurs for a number of 
quantum models [12]. From equations ()2.4[) and ()2.7|) . 

-^12 'S'12 Ti S12 T2 = S12 T2 3^2 Ti R12 ■ (2-15) 

Assuming the particular relative value of Ai and A2 such that R12 is proportional to a rank-one 
projector, the quantum determinant [H], [T3], is denoted by A, where 

R^2^ = RuSuTiS+2T2. (2.16) 

From equations ()2.1aj) and ()2.7j) 

Rl2^12 -^n,l Ln,2 = 'Jl2-^n,2 -^n,l -Rl2 ■ (2-17) 

Maintaining the same particular relative value of Ai and A2, the local quantum determinant is 
denoted by Qet(L„), where 

RuQ^^i^n) = R12S12 Ln,l Ln,2- (2-18) 

In [Appendix B[ it is shown that the quantum determinant factorizes in terms of the local quantum 
determinants as 

p 
A=n Qet(^n). (2.19) 

n=l 

In section [21 the quantum mappings of KdV type are considered. The quantum determinant and 
the local quantum determinants are central elements of the algebra for this model. Indeed, the 
quantum determinant will play a central role, in both the mathematical and conventional English 
sense, throughout the rest of this paper. 
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3. Quantum mappings of KdV type 

In operator form the quantum mappings of KdV type, as introduced in jSHj and |2I], are 

V2J-1 = V2j V2j = V2J+1 - aV2j~'^ + aV2j+2~^ (j = l,---,P), (3.1) 

with imposed periodicity condition fi+2P = Vi, P E N. The dynamical variables, Vn, are 
Hermitian operators, a is a real number parameter. In the notation of [201 the commutation 
relations of the dynamical variables read 

[vj , Vf] = h {5jj^+i - 6^+ij>) , (3.2) 

(h = —ih, where h is Planck's constant divided by 27r). The periodic initial value problem, from 
which the mapping arose, imposes on the mapping ()3.1|) the Casimirs 
p p 

5Z^2i = $^^2i-i =: i^ , (3.3) 

i=i j=i 

in such a way as to leave the value of these Casimir operators as a free parameter (it can easily be 
seen from the commutation relation, equation ()3.2j) . that this is a Casimir) hence, in the classical 
limit, we obtain what could be called a (P — l)-dimensional configuration space generalization 
of the McMillan mapping |T^. We assume that z/ 7^ 0. 

We need to point out that, at this stage, we are only concerned with the algebraic structures 
behind the integrability of the quantum discrete-time systems. Hence, as far as this paper is 
concerned, we will deal with operators, such as the {vk} := {vk}k=i...2P, on a strictly formal level 
(in the spirit of related work ^Oj)- This involves, for instance, assumptions on the invertibility of 
the operators, disregarding, for the time being, questions concerning the domains of the Hilbert 
spaces on which they act (we aim to return to this latter issue in subsequent publications). 

The Lax description of the mappings of KdV type is as follows: 

L, = V2, V2,^i V, = (^ l) (3.4) 

and X2J = A, X2J+1 = A + a. The associated monodromy matrix, T(A), is obtained by gluing the 
elementary translation matrices Lj along a line connecting the sites 1 and P + 1 over one period 
P, namely 



^(^) - [ C\7) o\7, ) - n ^»(A) ■ (3.5) 

The monodromy matrix has a natural grading in terms of the spectral parameter. A, 



A^ + A^-Mp_i + ... + A0 A^-iPp_i + A^-2Pp_2 + . . . + Pi 
X^Cp + X^-^Cp-i + ... + XC1 A^ + X^-^Dp^i + ... + XDi 







Tw = .p' ,\p-r^ ;■■ , " .p:\p-r. \'' ;v. ' • (s-e) 



Observe that A{X) and -D(A) are both monic polynomials in A. The time evolution is given by 
T(A) = M(A)T(A)M(A)-i M„ = | "^^ J | , (3.7) 
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where M is Mi, the discrete-time evolution operator at lattice site 1. More explicitly, bearing in 
mind that we are dealing with noncommuting operators, this gives us, 

wB + D J {wA + C - wBw - Dw) 
XB A-Bw 



T(A) 



where w := wi. As well as the mapping ()3.1|) . the Zakharov-Shabat condition 

a 



w„ 



V2n-\ ^ • 

V2n 



(3.8) 

reveals that 
(3.9) 



For the mappings of KdV type, the realization of the R and S matrices, which are solutions 
of the compatibility relations ()2.6aj) and ()2.6&|) under the condition ()2.7|) . is as follows: 



R 
R 

S- 



12 



12 



-"^12 "^12 



S- 



12 
Pl2 



12 



h 

a; 



Ai — A2 

F®E 



(3.10) 



S 



12 



s. 



21 5 



where the permutation operator P12 and the matrices E and F are given by 

/ 1 \ 



12 



10 
10 
1 



E 





(3.11) 



The realization ()3.10|) is assumed throughout the rest of the paper. 

Classically equation ()3.7j) gives us that the trace of the monodromy matrix is invariant under 
the discrete-time evolution. This argument no longer holds in the quantum case, as some of the 
matrix entries consist of noncommuting operators. As stated in the previous section, within the 
quantum case the invariants are given by expanding equation ()2.1H) in powers of the spectral 
parameter, A. For mappings of KdV type equation ()2.12|) gives 



K{\) 





1 + f 



thus 



r(A) = A{\) + 1 + - D{\). 



(3.12) 



(3.13) 



Equation (j3.10|) shows that R12 
A2 = Ai + h. Therefore, with Ai = 



is the fully antisymmetric projector, 1 — P12, when 



A 



A- and A2 



A + 



gives the quantum determinant for this model. It may be expressed as 



A(A) 



A_ 



D(A_)A(A+)-i?(A_)C(A+). 



A+, equation ()2.16|) 



(3.14) 



There are, of course, other equivalent expressions which can be obtained using the algebra (j2.4|) . 
for instance 



A(A) = ^A{K)D{\_) - ^B{K)C{\.). 



(3.15) 
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then it can easily be shown using ()2.18|) that the quantum determinant of the algebra ()2.1aj) can 
be written as, 

Qet(L„(A)) = ^d„(A_)a„(A+)-6„(A_)c„(A+). (3.17) 

Hence, from ()3.4p . Qet(L„(A)) = A+(A_|_ + a) and, therefore, from ()2.19p . 

A(A) = A^(A+ + a)^, (3.18) 

which manifestly belongs to the centre of the algebra. 

4. Sklyanin variables 

Following Sklyanin [2^] the operator zeros of B{X), x„, provide the separated canonical variables. 
By "operator zeros" it is meant that 

p-i 
B{X) = Bp ^11 {X-x^) , (4.1) 

71=1 



where the {xj} := {a;j}j=i,..p_i mutually commute. (In [Appendix C it is shown that Bp_i is 



equal to the Casimir ()3.3|) . We assume the mutual commutativity of the {xj}, this is consistent 
with equation ()2.4|1 .) Conjugate variables to the Xn are introduced by making the definitions 

X- = A{X) X+ = f 1 + y) DiX) (4.2) 

where the operator ordering prescription throughout this paper is that Xn is substituted for the 
spectral parameter. A, from the left, thus, 

X~ = Ao + XnAi + ... + x^, (4.3a) 

X^ = hDi+Xn{Di + hD2) + ... + x^. (4.36) 

As in [2^ the full set of commutation relations between these operators follows from the Yang- 
Baxter structure, equation ()2.4|) . and reads, 

[xm, Xn] = Wm,n (4.4a) 

X^Xn = {xn±h Smn)X^ M m,n {4: Ah) 

[X^,X^]=Q \/m,n (4.4c) 

[X-,X+]=0 Wm^n (4.4c?) 

X^X^ = a(x„±^) Vn (4.4e) 
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where A(A) is the quantum determinant of the model, given exphcitly in equation (j3.18p . We 
also have 

r(x„)=X+ + X-, (4.4/) 

which leads to the linear finite-difference spectral problem known as Baxter's equation (see, for 
example, ^ and |3Qj). 

The derivation of the Sklyanin algebra relations (equations fj4.4ffl|) to ( |4.4/| )) proceeds in the 
same way as in jSHI, but is slightly more involved as the initial equations from the Yang-Baxter 
equation are more complicated (there are extra terms due, essentially, to the non-ultralocal 
nature of this algebra). Remarkably, as the derivation is carried out, the extra terms vanish, 
leaving the Sklyanin algebra relations. 

The proof of the preservation of equation ()2.4|) . for both the temporal and spatial evolutions, 
is given in [Appendix A As a consequence we have the following result for the Sklyanin algebra 



under these evolutions. 

Proposition. The Sklyanin algebra relations are preserved under both the temporal and spatial 
discrete evolutions. 

The (extended) Yang-Baxter structure of section |2l leads efficiently to the preservation of 
the Sklyanin algebra relations, as is expressed in the proposition. The Sklyanin algebra relations 
are the real starting point of this work. We now turn to the equations of motion in terms of the 
Sklyanin algebra variables. 

5. Quantum discrete Dubrovin equations 

The equations of motion, for the temporal and the spatial evolutions, are derived in this section. 
The aim is to establish these discrete evolutions in terms of the Sklyanin algebra variables; this 
requires the reconstruction of w and one of the original dynamical variables, 1/^2, in terms of 
the Sklyanin algebra variables. Hence the issue of the reconstruction of the original dynamical 
variables in terms of the Sklyanin algebra variables is also, necessarily, addressed. 

In section 15.11 the invariants of both the temporal and spatial evolutions are expressed in 
terms of the Sklyanin algebra variables. 

In section 15.21 the reconstruction of w in terms of the Sklyanin algebra variables is given. 
This paves the way for the derivation of equations ()5.14j) and ()5.2()|1 . As is illustrated for the 
P = 2 and P = 3 cases in section IHl these two equations, along with the preservation of the 
invariants, give a well-defined temporal evolution. Hence, equations ()5.14|1 and ()5.2()|1 of section 
15.21 are what we mean by the temporal quantum discrete Dubrovin equations. 

In section EISl the reconstruction of l/f2 in terms of the Sklyanin algebra variables is given. 
Along with the reconstruction of w this paves the way for the derivation of equations ()5.29|) and 
()5.31|) . In an exactly analogous fashion to ()5.14p and ()5.20|) in the temporal case, equations ()5.29|) 
and ()5.31|1 give a spatial evolution. Hence, equations ()5.29|1 and 1)5.311) of section 15.31 are what 
we mean by the spatial quantum discrete Dubrovin equations. We conjecture that the spatial 
evolution allows for a reconstruction of all of the original dynamical variables in terms of the 
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unshifted Sklyanin variables. This is illustrated in the P = 3 case in section IHl and it is not 
technically difficult to con&m this for the next few larger-period cases. However, the calculations 
quickly become very cumbersome as the period increases. 



5.1. Invariants 

In this section expressions are given for the invariants, which are the coefficients of the various 
powers of A in equation ()3.13|) . in terms of the Sklyanin algebra variables. From the form of the 
invariant in terms of entries of the monodromy matrix, ()3.13|1 . and their gradation ()3.(i|l . 



TiXr. 



„P-1 



„P-2 



2< + <~ Ip-1 + <~ Ip-2 + .-.+ Xnh + h- 



(5.i: 



It is easily shown that /p_i is a Casimir operator. Taking into account the gradation given in 
equation ()3.6|) . equation ()3.14|) gives us that the quantum determinant 

A(A + |) = A^^ + A'^-^(P/i + Ap_i + Dp_i + h- Bp^iCp) + .... 
However, for the mappings of KdV type we have (from equation ()3.18|l ). 



h. 



A(A + -) = A^^ + y'^-\2Ph + Pa) + . . . . 



Therefore, using that the central elements Pp-i = Cp = v (as derived in Appendix C ) 

/p_i := Ap_i + P)p_i + /i = i^2 + P(a + /i). 
Now observe that, 



(5.2) 



h 



/I 
1 



X\ 
X2 



P-2 

P-2 
X2 



\ 



( t{x\) — 2xf 
p 



„p-i 



r\X2) 



x\ ^Ip-x 



1 
2x^ — X2 ip-i 



^X p 1 



„p-l 



.(5.3) 



X J 



iip-1 ) 



\ Jp_2 / \ 1 xp_i . . . Xp_i / \ r(xp_i) 

All terms on the right hand side are known (remember that t(x„) = X+ + X^). Equation ()5.H) 
may also be rewritten as 

-1 



h 
\ Ip-i J 



( x\ 

X2 



XI 
2 



P-l 



\ 



P-l 



„P-l 



( 



T\Xi] 
t{x2) 



2xf - Jo 
2x1 - h 



\ 



\ ^(a;p-i) 



-^x p_1 



(5.4) 



hi 



\ Xp_i X\>_^ . . . Xp_^ I 

Note that within the VanderMonde matrices all of the entries commute. 

As the quantum discrete Dubrovin equations, which are derived in the following sections, are 
also matricial equations, it is expedient to introduce some specialized notation. The symbol M. 
is introduced to denote the VanderMonde matrix, the symbol T> to denote the diagonal matrix. 



M 



( Xi 

X2 



X, 



Jbn 



X 



p-l 



\ 



Jjg 



„P-i 



\ a;p_i Xp_^ 



V 



Xp_i ) 



/ xi 

X2 

V 



\ 


Xp_i j 
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Vectors are indicated by a bold typeface, and all are P — 1 dimensional. The vector e = 
(1, 1, ... , 1)*. The vector consisting of ordered entries, labelled n to P — 2 + n, for any integer 
n, is denoted by a bold typeface with a subscript n, for instance yo = {yo,yi, . . . ,yp-2y, 
yi = (l/i,Z/2, 



. . ,|/p_i)*. In this notation equations ()5.3p and ()5.4|) are rewritten as 
lo = M-^V (X+ + X^- - V^2e - V^-^Ip^^e) 



and 



Ii = M-^ (X+ + Xf - V^2e - Joe) , 
respectively, where Iq and Ii denote the left-hand sides of ()5.3p and ()5.4|) respectively. 



5.2. Temporal equations 



In section 15.2.11 we give the reconstruction of w in terms of the Sklyanin algebra variables. For 
heuristic reasons the matricial equations are written out in full in this section. In section 15.2.21 
the temporal quantum discrete Dubrovin equations are derived. 



5.2.1. Reconstruction of w. Observe that, from the definitions ()4.2j) (and remembering that Xn 
are substituted for A from the left), 

:= D{xn) = -""'^ 

A — Xji. 



D{\) 



-X„ 



Xn + h 

(We take the casual attitude of writing 1/Y for the quantum inverse of the operator Y .) Hence, 
still working from the definitions. 



/ ^1 \ 

A2 

V Ap-1 1 

f D, \ 

D2 



\ Dp-i I 
So, it is seen that. 



/ 



X\ 
X2 



Xi 


X? . 


.. xr 


X2 


iX'o 


. . %JjO 



1 



X 



p~l 



1 



djr) 



\ Xp-i Xp_i 



\ 



.P-l 



( ^r 



^p-i / 



<-' 



XT, 



V^P-i 



p 

a;f 

p 

x^ 



^0 



\ 



X 



p-\ 



Ao J 



X. 



p-1 



„p-l 



-1 / ^^X+ - T^ 



Xp_i I 



h+X2 2 



Xn 



\ h+xp- 



^^p-l -^P-l I 



(5.5) 



(5.6) 



/ 



A2-D2 



hD2 



\ 



\ Ap^i - Dp_i + hDp ) 



( x\ 

X2 



SjA 



JjO 



\ ^P-\ Xp_^ 



Xi 



n 



-1 



„p-l 



Jun 



X 



p-i 



( 



\) 



X- 



X 



p-1 



h—x 



^X+ - Jo \ 



h+. 
h—X2 y+ 



h+x: 



x^ 



h+xp_i^P-l 



'"J 



(5,7) 
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The time-update of the inonodromy matrix, equation (j3.8|) . and a commutation relation contained 
within ()2.10|) give 



A = Bw--{A- Bw)+ (l + -] D 



D = A-Bw. 



Substituting a;„ for A from the left, and again using the definitions (j4.1|) and (|4.2p . gives, 

XnDi + xl^D2 + . . . + X^ = X~. 

From here we easily obtain (using Aq + hDi = Iq = Iq) that 



(5.8) 

(5.9a) 
(5.96) 



/ 4i - ^1 + hD2 \ 
A2-D2 + hD-i 



( a;i 
X2 



X -| 



Jun 



JL 1 



„p-i 



\ 



-1 



X 



p-l 



\ Xp_i Xp_^ 



( ^i+-^r 



^p-i / 



-^2 ~ ^2 ~ -^0 



\ ^p-i - ^p-\ - h J 



(5.10) 



\ v4p_i - L)p_i + /iDp / 
Equation (j3.8|) also gives 

1(A) - D{\) + A{X) - D{X) = wB{X) + B{X)w. 
A consideration of the leading term as A — > 00 shows that 

Ip_i - 5p_i + Ap_i - Dp_i = 2Bp^iw, (5.11) 

as Bp_i is actually the Casimir (|3.3p (this is shown in Appendix C ). Therefore 

w = — (Ip-i - 5p_i + Ap_i - Dp_i). 
2v 

Expressions for Ap_i — -Dp_i and v4p_i — -Dp_i in terms of the Sklyanin algebra variables, 
{xi,Xf} := {xi,Xf}i=i p_i, follow from equations ()5.7p and ()5.10|) . To obtain an expression 
explicitly in terms of {xi.Xf}, Iq may be replaced by using (j5.4j) : then Cramer's rule, along 
with ()5.2j) for the value of the central element /p_i, allows one to deduce that, in terms of the 
Sklyanin algebra. 



w 



V 



p-\ /P-i 



.^+pa+(p-i)/.+(-i)-$:(n-4 






Xi Xr. 



2x 



p-i 



h + Xn " 



-X„ 



Xr. 



.(5.12) 
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5.2.2. Temporal quantum discrete Dubrovin equations. Equations (j5.7p and (j5.1(jp give, in the 
notation introduced in section I^TTl 

M-^ fx^ + {hi + V)-\hl - P)X+ - Joe) = M-^ fx+ - X^ - Joe) . (5.13) 



This constitutes the (temporal part of the) quantum discrete Dubrovin equations, as given 
classically in ^^I- Note that the time-evolved variables still obey the relations ()4.4ajl to ( 4.4j| ) 
(as stated in the proposition of section Hj). Along with the time-update invariance of ()5.4|1 . ()5.13j) 
leads to 

M-^ (v^e - {hi + V)-^VX+) = M-^ (v^e - X^) . (5.14) 

The equations for the elementary symmetric polynomials in {xi} in terms of {xj, X^ } will 
now be obtained. A consideration of equation ()4.H) shows that these follow immediately from 
the coefficients of different powers of A in B{\), that is, from {Bq, . . . , Bp^2}- Equation ()3.8|) 
and a commutation relation contained within ()2.1()|1 give 



B(X) = - 



Aw - Bw' +[1 + -]{C- Dw) 



(5.15) 



Hence, to obtain the elementary symmetric polynomials in {xj} in terms of {xi,Xj^}, C{\) 
must be expressed in terms of {xi^X^}. From equations ()3.14|) and (j3.18|) for the quantum 
determinant, 



B{\-h)C{\) 



X 



X-h 
Consider the first term on the right-hand side 

A 



D{X-h)A{X)-X^{X + af. 



2P-1 



X-h 



D{X - h)A{X) = ^ A^+^ J2 



j=Q 



ft=0 
k>j-P 
k<P~l 



P-l 

E 

J=k 



[-hy-'Di+, 



A,. 



(5.16) 



(5.17) 



where Dp = 1 and (J denotes the binomial coefficient, (^) = /!/[(/ — A;)! A;!]. Whence, substituting 
x„ for A from the left, 



C [Xr, 



Bp-i 



p-i 

n 



J-'r. 



1=1 

2P-1 



Xi — h 

E-;^' E [E 

j=0 



fe=0 

k>j-P 

k<P-l 



{-hy-'Di^, 



Aj~k ~^ ^n [-^n + 0,) 



(5.18) 



A consideration of the A° term of D{X) in equation (j5.8|) reveals that Aq = Bqw. As w is given in 
terms of the Sklyanin algebra, {xi,Xf}, in equation ()5.12|) . we have, along with equations ()5.5|) 
and (j5.6p . all of {A„} and {Dn} in terms of {xi,X^}. Therefore equation (j5.18p gives C{xn) 
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entirely in terms of the Sklyanin algebra. Consider now equation (j5.15|) with x„ substituted for 
A from the left, this defines, 






(X--X^)w+[l + -]CiXr. 



(5.19) 



The right-hand side of ()5.19p may be expressed entirely in terms of {xi,X^ }, for brevity 

denote the right-hand side, strictly in terms of {xi,Xf}, by i3„. Therefore, with Si = 
{Bi, B2, . . . , Bp^iY, we have the following expression. 

Bo = M-^V (Bi - V^-^Bp^ie) , (5.20) 



which gives the elementary symmetric polynomials in {xi} in terms of {xi,X- }. 

Equations (j5.14p (or (j5.13|) ) and (J5.2(J|) constitute the temporal quantum discrete Dubrovin 
equations. 

5.3. Spatial equations 

In this section we consider an evolution along the lattice at a constant time level. Specifically, 

T = LpLp„iLp_2 . . . L2L1 h-i> T = LiLpLp^i . . . L3L2 !—>■■■ (5-21) 

(The hat, ^ , is used to denote the evolution along the lattice at a constant time-level.) In 
the temporal case the reconstruction of w is required for the quantum discrete Dubrovin 
equations. The spatial equations also require the reconstruction of 1/^2, this is achieved in 
section 15.3.11 Together with equation ()3.9p . this also gives the reconstruction of Vi. If one 
requires a reconstruction of the original dynamical variables, {ffc}, in terms of the Sklyanin 
algebra set up at each dynamical variable's particular lattice site (that is, Vi and V2 in terms of 
{xi,Xj^}, f3 and v^ in terms of {xi^Xf"}, etc.) then the reconstruction is complete. However, 
we conjecture that the spatial quantum discrete Dubrovin equations, as derived in section I5.3.2[ 
allow for a reconstruction of all of the original dynamical variables in terms of the unshifted 
Sklyanin variables {xi^X"^}. 



5.3.1. Reconstruction of ^. In [Appendix C it is shown that -Bp-i is equal to the Casimir ()3.3|) . 



it is obvious that -Bp-i also is (as it is still the sum over all {f2j})- Now, from the definitions of 
the "conjugate variables" in section EJ 

Ai - Di + KD2 ± 2Bo— = M-^ f X7 + ihl + VV^ihl - r')X+ - Joe t V^2Bp^i—e\ .(5.22) 

Spatially updated terms follow from ()5.2H) (or, equivalently, ()2.14|) ). Using the commutation 
relations in (j2.13p . 

A + ^D = -^{A - Bv{){v2Vi + x + a)+(l + ^ T^(^ " ^^1)^2, (5.23a) 

A A + a \A/A + a 

D = -^[-Av2 + B{\ + a + viV2)]vi + ( 1 + ^ ) T^[-Cv2 + D{\ + a + V1V2)], (5.236) 

A -\- a \ X I A -\- a 
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B 



X{X + a) 



[-Av2 + B{X + a + t;it;2)](wiW2 + h + X + a) 



h 

+ 1 + T 



So, 



^ h^ - - 1 

A+-D-D + 2XB— 

A V2 



Xj X{X + a) 

X + a + 2h 



Cv2 + D{X + a + f it;2)]'i^2- 



(5.23c) 



A+[1 + ^]d + 2^^±^±^i?(A + a + v^V2)-. (5.24) 
A + a \Ay X + a ^ V2 



From equation ()5.24j) one obtains 



Ai - Di + /iDs + 2Bo — 

V2 



= M-^ ( X+ - (D + al)-\V + al + 2hl)X]; - IqB - D^25p_i— e J 
It is easily shown that the leading term in equation ()5.24p as A — ^ oo gives 



(5.25) 



2 (bp_,w + fip_i- ) = (lp_i - 5p_i + h + 2Bp_2—) + (Ap^i - Dp_i + h- 25p_2-).(5.26) 

V ^2/ V2 V2 

Equations (J5.22J) and (J5.25J) along with (j5.7p . (j5.1(Jj) . and (j5.11|) from the reconstruction of w, 
and Cramer's rule, lead to 



1 _ 1 

V2 V 



/p-1 



P-1 . /P-1 . 






\ 



/ 



x„ + a 



-x: 



(5.27) 



Along with w, as reconstructed in equation ()5.12p . and equation ()3.9p for Wn in terms of the 
dynamical variables {ffc}, equation ()5.27p leads immediately to the reconstruction of v\. 

5.3.2. Spatial quantum discrete Dubrovin equations. Equations (|5.22p and (|5.25|) give 



M-^ fxr + [hi + V)-\hl - P)X+ - Joe - V^2Bp_i — 

\ V2 

M'^ fx+ -(V + al)-\V + al + 2hl)X7 - Joe - V^2Bp^i—e) , (5.28) 



for the spatial part of the quantum discrete Dubrovin equations. As the invariants of the time 
evolution are also invariants of the spatial evolution we may use them, in the form given in 
equation ()5.4|) . to rewrite ()5.28|) as 



M-^ {hi + Vy^VXf - P^ 1 - Bp_i— e 



V2 



M-^ ({V + al)-\V + al + /il)X^ - V^ (l - 5p_i— j e) . 



(5.29) 
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The expressions for the elementary symmetric polynomials in the "spatially updated" 
variables, {xj}, in terms of {xi,Xf} are now given. Upon substituting x„ for A from the left, 
equation ()5.23c|l becomes 

1 + -) , ^ ^ M xn)vl. (5.30) 

Note that C{xn) was expressed in terms of {xi^Xf} by equations ()5.5|1 . ()5.(jp . and ()5.18j) : w was 
reconstructed in terms of {xj, Xf} in equation ()5.12|) : t>2 was reconstructed in terms of {xi, Xf} 
in equation ()5.27j) . Therefore the right-hand side of ()5.3()j) can be given entirely in terms of 
{xi, Xf}. For brevity denote the right-hand side of ()5.3U|) expressed strictly in terms of {xj, Xf} 
by Bn- With this notation 



Bo = M-'V ( Bi - V^-'Bp^ie] , (5.31) 

giving the elementary symmetric polynomials in {xi} in terms of {xi,Xf}. 

Equations ()5.29p (or ()5.28|) ) and ()5.3H) constitute the spatial quantum discrete Dubrovin 
equations. 

6. Examples 

The well-defined temporal evolution which follows from the quantum discrete Dubrovin equations 
is illustrated for the P = 2 case in section 16.11 and for the P = 3 case in section 16.21 
The reconstruction of the original dynamical variables, {vk}, in terms of the Sklyanin algebra 
variables, {xi,X- }, is also performed in the P = 2 and P = 3 cases. The reconstruction is 
achieved, essentially, via a "spatial evolution" using the spatial part of the quantum discrete 
Dubrovin equations. 

6.1. The P = 2 case 

6.1.1. Temporal evolution With P = 2, equation ()5.12j) gives 

w = -(u^ + 2a + h + 2x ^--X+ - -xA . (6.1) 

u \ X + h X J 

Equation ()5.14|1 gives directly that 

X ^X+ = x--X-. (6.2) 

X + h X 

Along with ()6.1|) . equation ()5.20|) leads to 

i = i„f,_ix-V (6.3) 

z/ V X 
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It is easily seen that 

1 + - ) -^X+ = --X- + -lo + /i + 2x. 
X J X + h X X 

Hence we may consider the evolution given by the quantum discrete Dubrovin equations to be 
that of X and ■^^X'^ (or ^X") along with the preservation of the invariant, Iq = Iq . . .. Equation 
(JE21) gives Ji = 1/2 + 2(a + h). 

6.1.2. Reconstruction Setting P = 2 in equation ()5.27j) gives 

- = i (l - ^^X-) . (6.4, 

V2 v \ x{x + a) J 

It follows trivially from the Casimir (j3.3|) . f 2 + f 4 = v, that we also have v^ in terms of (x, X^). 
Equation ()3.9j) along with (jfj.lj) and ()(i.4|l gives the reconstruction of fi, and the Casimir ()3.3|1 . 
"^1 + "^3 = 1^1 gives the reconstruction of ^3. Therefore, full reconstruction has been achieved. 

6.2. The P = 3 case 

6.2.1. Temporal evolution For this section we define 



2^2 — a^i V^i + ^ X2 + h J X2 — Xi \Xi + h X2 + h 



a^2 — a^i \Xi X2 J X2—x\ \x\ X2 

Then, with P = 3, equation ()5.12|) gives 

w = - [z/2 + 3a + 2/i + 2(xi + X2) + VTa + ^2] • (6.5) 

Equation ()5.14|) gives directly that 

X1X2 + TTi = a;iX2 + 01, (6.6) 

^1 + ^2 + 7r2 = a^i + a;2 + 02- (6.7) 

Along with ()6.5|) . equation ()5.20|) leads to 

x\X2 = -w (xia;2 + 0i) , (6.8) 

5i + £2 = -ui (xi + X2 + 02) ^[2xiX2 + h(xi + X2) + h"^ - 3a2 

- (Xi + X2 + IK) (Xi + X2 + 2/1 + 3a) + VTi + 01 + 7r202] . (6.9) 

It is easily seen that 

h 1 

-TTi + 7r2 = -02 H h - h- 2(xi + X2) 



X1X2 X1X2 
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and 



1 + 



h{xi + X2) 
X1X2 



TTi - hTT2 = -01 H Iq + h - 2X1X2. 

X1X2 



Hence we may consider the evolution given by the quantum discrete Dubrovin equations to be 
that of the elementary symmetric polynomials X1X2 and Xi + X2, and tti and 7r2 (or 0i and 02), 
along with the preservation of the invariants, Jj = Jj . . .. The form of the invariants follows from 
(Q and (Q. 



6.2.2. Reconstruction Setting P = 3 in equation ()5.27|) gives 



1 _ 1 

V2 V 



1 + 



-^1 -, r^2 

X2(X2 + a) 



(6.10) 



In the P = 3 case, equation (j5.29p of the spatial part of the quantum discrete Dubrovin equations 
reads 



-X1X2 



\^ V2 J X2—X1 yxi+h i X2+n ^ 

(xi + X2) (1-^) + ^-^-^ (^X+ - ^X^ 

\ i- ' ' \ ""2 J X2-X1 yxi+h i X2+h ^ 

-X1X2 f 1 - -) ^— [^ f 1 + -^) X7 -'^{1 + -\-\ X, 

\ 112 J X2—x\ yx\ y x\+a J i X2 \ X2+a J ^ 

(xi + X2 - /i) f 1 - -) + -^— {-X7 - ^Xo 



The top — (a + /i)bottom of ()6.11|) leads to 



. (6.11) 



[-X1X2 - (a + /i)(xi + X2)] ( 1 \~- \ - , , 

V2 X2 - Xi V Xi + /l 



X2 + a + /i ;r>+ xi + a + /i -^ 

^1 ^ ; — ^2 

X2 + /i 



a^hy\\--\. 



V2 



Therefore, 



1 _ 1 



1 + 



:^t- 



(6.12) 



-X^ 



3^2 — 3^1 V (^1 + /i) (a^i + a + /i) ^ (a^2 + ^) (3^2 + a + /i) 
It follows trivially from the Casimir (j3.3|) . f + f 2 + ^^4 = '^j that we also have ^4 in terms of 

{Xj,Xj }j=l,2. 

The reconstruction of {f2j+i} follows from that of w. It follows from equation ()3.9j) that 
v\, in terms of {xi,Xj^}j=i^2, is obtained from equations ()6.5|) and ()6.10|) . Using the form of Jq 
which follows from ()5.3p . one obtains. 



w 



V 



Jo -h TTi 



X1X2 



X1X2 
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Now, as Jo is invariant under spatial updates, it may be written, using the same expression from 

equation ()5.3j) . at any spatial level. In simple terms, as Jq = /q = /o = • • •, the same expression 
for Iq may be used with any number of hats above or below the operators. Hence, 

h ^ 



w 



1 



-In-h 



-TTi 



X1X2 a;iX2 

The top equation of the spatial quantum discrete Dubrovin equations, (jfj.llj) . gives 



(6.13) 



xia;2 



-TTi 



X1X2 



a;ia;2 
1 



V 

1 

^2 



X\X2 X2 — Xi 

and equation ()5.31|) gives. 



^11 

Xi 



h 



Xi + a 



AT- 



-(1+ 

X2 V ^2 



h 



Xo 



(6.14) 



X1X2 



(XiX2 + 0l) V2W + -X1X2V2W H (XiX2 + TTi) ■U2 

az/ a az/ 

h 



au^ 



-\^x\X2 + /i(xi + X2) + /i — 3a 



(xi + X2 + 2h){xi + X2 + 2h + 3a) + vti + 0i + ■rc2(j)2]v2- 



(6.15) 



So, with equations ()6.5|) and ()6.10|) . this gives £1X2 in terms of {xj,Xj^}. Therefore we 
have reconstructed w in terms of {xi, X^ }. With the reconstruction of V4 this then gives a 
reconstruction of V3, and from the Casimir ()3.3p . f 1 + ^3 + f 5 = z^, this gives a reconstruction of 
f5. Therefore, full reconstruction has been achieved. 



7. Conclusion 

The quantum discrete Dubrovin equations have been derived. Equations ()5.14|1 and (j5.2()j] give 
a temporal evolution and ()5.29|) and ()5.31|) give a spatial evolution. 

The classical discrete Dubrovin equations were published in [iJI and J22], but in jT^j the 
focus was on temporal equations only. In the temporal case, the work presented here is very much 
analogous to the classical case. The classical {h — * 0) limit of ()5.13|) is the discrete Dubrovin 
equations as presented in j^ . This is seen as follows. The derivation of the (classical) discrete 
Dubrovin equations employs the invariant spectral curve 

det(T(A) -v)=0, (7.1) 

which defines a hyperelliptic curve of genus g = P — 1. The classical equations are written 
in terms of the discriminant of the hyperelliptic curve ()7.1|) . R{X). This discriminant may be 
expressed as 

1 



i?(A) = -(A(A) + D(A))^-det(T(A)). 



(7.2) 
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Following [T7j, we see that (j7.2|) implies that 

1(A(A) - D(A)) = ^^R{X)-BiX)CiX), 

where k denotes the sign ± and corresponds to the choice of sheet of the Riemann surface, the 
condition k = —k, can be seen by also considering ()H.8|1 and evaluating at the (no longer operator) 
roots of -B(A). A quantum deformation of the expression ()7.2|) . evaluated at an operator root of 
B{X), Xn, is 

^(^n) = \r{xnr - I (A(^n ' ^) + A(x„ + ^^j = \ {X' " X^)' ■ (7-3) 

Therefore, in the classical limit, the quantum discrete Dubrovin equations, ()5.13|) . become 

M-' (^k^/K(^ - \he) =M-'(^- K v/R(^ - ^/oe) . (7.4) 

where a/R(xi) = {^^/R(xl), . . . , a/-R(xp_i))*. Equation ()7.4j) is the classical discrete Dubrovin 
equation, as first presented in [T7|. (The notation of [T7j is such that the roots of B{X) are 
denoted by {fii}, rather than by {xj}.) Classically, the discrete Dubrovin equations are the 
intermediate step towards the parameterization of the orbits of the classical map. In ^j the 
parametrization of the solutions of ()7.4j) in terms of Abelian functions of Kleinian type was 
discussed, and illustrated in the P = 2 and P = 3 cases. From the new perspective of this paper 
it is surprising that the classical limit of ()5.20|) is not required for this parametrization, and, 
indeed, does not feature in [T2j. 

Returning to the quantum equations, the next issues to be addressed concern the 
representation theory. Then, drawing inspiration from the classical discrete Dubrovin equations, 
we would hope to be able to construct explicit expressions for the quantum propagators 
interpolating over an arbitrary number of discrete-time steps. Effective mechanisms for 
computing expectation values and long-time asymptotics for the transition amphtudes would 
also be desirable corollaries of this proposed work. Extensions of the present work to the higher 
rank case associated with quantum mappings in the Gel'fand-Dikii hierarcy ^Sj, ^H], also form 
the subject of future work. 
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Appendix A. Quantum invariants 

Consider the following evolution of the monodromy matrix: 

r = NTN^^ (A.l) 

(so, to be specific, N could be Mi or Li). In this appendix it will be proven that if N satisfies 
the relations 

i?+ N1N2 = ATaATi/?- (A.2) 
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and 

TiATfi^+ATa = N2S^^TiN{^ (A.3) 

then the Yang-Baxter relation p.4|l is preserved through this evolution and, moreover, 

r(A) = tr(K(A)T(A)) 

is invariant under this evolution, for a certain K{X) which will also be derived. First the 
preservation of the Yang-Baxter relation ()2.4|1 under the evolution (jA.lll is shown. It requires 
only ()A.2|1 . ()A.3|1 . and the constraints on the R and S matrices that -R^^'S'il = >S'j^-R5''2 5 

Rt,T[ St^n = Rt2Ni T, N^' St^N^ T^ N^' 
= Rf^Ni N2 S^^Ti N^' T2 N^' 
= N2 N, 5+ i?^2^i ^r' T2 N^' 
= N2 Ni St2Rt2Ti St2T2 N2' A^r'^iV' 
= N2 Ni S21T2 SfaTi i?r2^2"' ^T' V' 
= T^ St,N, N2 S^,T^N{' N^'St2"R^, 
= T^ S^^T[R-,, . (A.4) 

The derivation of the invariants of the evolution ()A.1|) will now be given. Introduce the tensor 
Ku = P12K1K2 , (A.5) 

where P12 is the permutation operator, which satisfies the relations 

Pi2{A®B) = {B®A)Pi2 Pi2 = P2i tr2Pi2 = li. (A.6) 



Choosing Ai = A2, we can take the trace of ()A.3|) contracted with K12. The left-hand side leads 
to 

tT^tT2{KuTiN{'St2N2) = tT2{K2T2N^'tT,{PuK2St2)N2) = ii{KT) (A.7) 

provided that 

tn{Pi2K2St2) = I2 . (A.8) 

Under the same condition, ()A.8|) . we have, from the right-hand side of equation ()A.3|) . that 

iiiii2{Kr2N2S^^TrN^^)= iir{K^N^ii2{Pi2KiS^2)TiN^^) = ii{KT') . (A.9) 

Hence, if K{\) is a solution of equation ()A.8|1 . we have the invariance of tr(K(A)T(A)) under the 
evolution given by equation ()A.1|) . 

A solution to equation ()A.8j) is found by taking 

K2 = tri{Pi2*^[(*^5+2)-']} . (A.IO) 

This is most easily verified by introducing the twisted product 

X12 * Y,2 :='' (*^Xi2*^Fi2) =*^ (*^ri2*^Xi2). 
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An inverse with respect to the * product is 

So, with equation ()A.10|) . 



tri(Pi2i^25+) = tritr3(Pi2P32(53+)-^*5+) 

= tri tr3(P32Pi3(S'32) *S^2) 

= tr3(P32 tri(Pi3(53+)-^*5+)) 

= tr3(P32 53+*(53+)-^*) 
= tr3(P32 I32) = l2- 

For the proof of the commutativity of the invariants r(A) = tT{K{X)T{X)) with the form of K{X) 
given in equation ()A.10|) we refer the reader to 



Appendix B. Quantum determinant factorization 

In the case of ultralocal models it is straightforward to show that the quantum determinant is 
equal to the product of the local quantum determinants of the constituent L operators. This 
factorization will now be proven for the present non-ultralocal case, which is less straightforward. 
Equations (ITTal) to (ETcl) and (|2IZ|) give 

mm mm 

-^12 11 Lj,l 11 Lk,2 =11 Lj^2 11 -^fc,l R12 5 (B.l) 

j=l k=l j=l fc=l 

for m < P — 1 . From equations ()B.H) and (|2.7p , 

mm mm 

-^i2'S'i2 II Lj^i II Lk^2 = S12YY Lj^2 II Lk^i R^2- (B-2) 

j=i fe=i j=i fc=i 

Therefore, for the particular relative value of the spectral parameters Ai and A2 such that R12 is 

proportional to a rank-one projector, equation ()B.2|) is defined to be equal to RuQ^^iYYjLi ^j)- 
Consider the left-hand side of equation ()B.2|) . 



m mm 

Pj"2Qet(|| Lj) = R^2S^2 II Lj,l II Lk,2 
j=l j=l k=l 



m— 1 m— 1 



Rl2^12^m,l Lm,2 ^21 || ^j,l || ^k,2 
j=l k=l 



m— 1 m—1 



^12 Lm,l Lm,2 Rl2^12 || -^j,l || ^k,2 
j=l k=l 
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m—l 

= QetiLjR^^Qetill L,,i), (B.3) 

and, therefore, 

n L, =n Qet(L,), (B.4) 

where m < P — 1. The proof of the same resuh via the right hand side of ()B.2j) follows similarly. 
The quantum determinant for the Yang-Baxter structure ()2.4|) . A, was introduced in 
equation ()2.1fij) . From the left-hand side of ()2.15|1 . 



p-i p-i 

= R12S12LP1LP2S21 J_J_ Lji J^J^ Lk2 

j=i k=i 

p-i p-i 

= S^2^P,2Lp,lRl2^12 II -^j,l II ^k,2, (B-5) 

j=l k=l 

(the same result follows similarly for the right-hand side) and hence, from equation ()B.4|) . we 
have 

p 

A=n Qet(Li). 
i=i 

Appendix C. -Bp 1 and Cp are in the centre of the algebra 

It can be shown, using the commutation relations ()2.4|) with the explicit realization of the R and 
S matrices ()3.1()j) and the gradation of the monodromy matrix ()3.fj|l . that -Bp_i and Cp belong to 
the centre of the algebra. The model specific proof given here, however, is more constructive for 
our purposes as it reveals -Bp_i and Cp to have the value prescribed to the Casimirs in equation 
()3.3|1 . Consider Lp+i, 

( \ + a + V2P+2V2P+1 V2P+2 \ /^ , X 

^""^ - 1^ Xv2P,, X ) ■ ^^-'^ 

From the definition of the monodromy matrix, equation ()3.5|) . and the grading of the monodromy 
matrix for period P, equation ()3.(ij) . along with equation (jdlj) . we observe that the monodromy 
matrix of period -P + 1, Tp^i, has the form 

A-^^^ + A'^[y4p_i + a + V2P+2V2P+1 + V2p+2Cp] + . . . X^[Bp^i + V2P+2] + ■ ■ ■ 

A^+i [Cp + V2P+1] + ... A^+i + A^[Dp^i + V2P+iBp_^] + . . 
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The B and C entries give us recursion relations for the coefficients of the highest order of A. Using 
these, in conjunction with the information from ()C.1|) that for the period = 1 case Cp = vi and 
Bp^i = V2, we obtain that for the period = P case 

p p 

Bp-l = ^t^2n Cp = "^Vin-l- (C.2) 

n=l 71=1 

For the mappings of KdV type these summations are the Casimirs given in equation ()3.3p . 
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